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Abstract 

Motivated by entropy estimation from chaotic time series, we provide a 
comprehensive analysis of hitting times of cylinder sets in the setting of Gibb- 
sian sources. We prove two strong approximation results from which we easily 
deduce pointwise convergence to entropy, lognormal fluctuations, precise large 
deviation estimates and an explicit formula for the hitting-time multifractal 
spectrum. It follows from our analysis that the hitting time of a n-cylinder 
fluctuates in the same way as the inverse measure of this n-cylinder at 'small 
scales', but in a different way at 'large scales'. In particular, the Renyi en- 
tropy differs from the hitting-time spectrum, contradicting a naive ansatz. 
This phenomenon was recently numerically observed for return times that are 
more difficult to handle theoretically. The results we obtain for return times, 
though less complete, improve the available ones. 
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1 Introduction 



The setting of this work is a model for the following situation: One has time series 
xi,X2,X3, ...,xn, ... assumed to be typical realizations of some dynamical system. 
We further assume that the invariant measure of the system is Gibbsian. We then 
want to evaluate the entropy, the Renyi spectrum and the hitting-timc/return-time 
spectrum of the system by using estimators based on hitting and return times. Our 
aim is to analyse the fluctuation properties of these estimators to have an a priori 
control on what we compute in practice. Regarding the estimator of the hitting- 
time/return-time spectrum, the issue is even to determine to what it converges. 

Before being more specific and describing these estimators, let us temporaly 
adopt a general point of view on hitting and return times. The first result in the 
early Ergodic Theory of dynamical systems is Poincare's recurrence Theorem, see 
e.g. or virtually any textbook on ergodic theory. Colloquially it states that in 
any dynamical system preserving a finite measure, typical orbits have a non-trivial 
recurrence behavior to each set of positive measure in that they come back infinitely 
often to it. More formally, let T : X ^ X he a, measurable transformation of the 
set X, the 'phase space', and /i a T-invariant probability measure on X. Poincare 
recurrence Theorem asserts that if C X is a measurable set of positive measure, 
then there is a subsequence {ni}i>i of the positive integers such that r"'a; € A for 
/i-almost every point x € A. One can also ask whether the orbit of a point a; ^ A 
will eventually enter A. A sufficient condition is ergodicity: for //-almost every x, 
there is a finite time hq = no(x) such that T""a; S A. This is the first hitting time. 

It is natural to seek for quantitative descriptions of this recurrence and hitting 
time. Kac's Lemma, see e.g. [21], states that the conditional expectation (first 
moment) of the first recurrence time to A is equal to l/fj,{A). This agrees with the 
intuition that the smaller the measure of A is, the longer it takes to come back. 
But ergodicity alone is not sufficient in general to guarantee the finiteness of the 
expectation of the first hitting time. Sufficient conditions on the mixing properties 
of the dynamics were recently given to ensure the finiteness of finitely many or all 
moments of hitting and return times [H] ■ We shall use them below. 

To go further, we turn our attention to sets A which are cylinder sets. Let A be 
a generating partition of X. To each point x G X we associate its n-cylinder, An{x), 
defined as the intersection of all the elements of A, T~^A, T'^'^^A containing 
X. Let us phrase the following remarkable result |33| : The return time r„(a;) of a 
/X- typical point x into n-cylinder An{x), behaves as follows^: 

r„(x) xexp(n/i^(r)) (1) 

where h^{T) is the measure-theoretic entropy of the system. An analog result is 
available for the first time w„(x, y) the orbit of x enters the n-cylinder An{y) about 
a point y, where x and y are 'randomly chosen' according to independently of one 
another: 

Wn{x,y) X exp{nh^{T)) (2) 

(But ergodicity is not enough: 'strong' mixing properties are necessary 1201 ■) 
mulas ^ and ^ give a simple entropy estimator based on the observation of a 
single typical orbit of the system. 

It is natural to look for the asymptotic law of return times, that is to study, for 
each t > 0, 

\ z e An{x) : rn{z) < (3) 

-"^The symbol X c,j' precisely moans lim,j^oo ^ logfon = lim„^oo ^ logc„ (/^ or ^x/x almost 
surely for pointwise quantities). 
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when n — !■ cx). An easy heuristic argument shows that one must rescale t by a 
quantity proportional to ^{An{x)) to obtain a non-trivial limiting law. Typically 

fi{Ar^{x)) ^ eM-nh^{T)) (4) 

by Shannon-McMillan-Breiman Theorem. 

Under an ever lessening set of hypotheses on the nature and the speed of mixing 
of the system, the limiting law has been proved to be the exponential law both for 
hitting and return times, see for instance |Hll8j. 

The purpose of this paper is to analyse the fluctuations of return times and 
hitting times by using approximations by the inverse measure of cylinders. Indeed, 
combining (Q) and or (0) and l@J we get 

Our aim is to sharpen these rough relations which are only pointwise. More pre- 
cisely, we want to compare the fluctuations of these quantities. Do they have the 
same log-normal fluctuations ? The same large deviations ? On another hand, there 
is a fundamental class of ergodic measures, namely Bowen-Gibbs measures [71 119|. 
such that ^ 

ji{Anix)) ^ e^p{Sn^p{x)) (6) 

where Sn^{x) is the 'energy' of the cylinder An{x). Therefore, one can reduce the 
study of logr„(a;) or logw„(x,?/) to that of 5„(/?(x), which is much more easy to 
tackle and, indeed, all is known on the fluctuations of Sn^{x). 

A related issue is to compute the so-called multifractal spectrum or Renyi en- 
tropy , defined as 

j fi{An{x))~'^ d^{x) X exp(nqM((;)) as n ^ oo . (7) 

In view of practical computation of the Renyi entropy, when one only has at hand 
a time series, it is tempting to make the ansatz 

^.{An{x)y^ ^ rn{x) 

in ^ and to evaluate the integral as a Birkhoff average. This was done by Grass- 
berger |16| . The implicit assumption is that these two quantities have the same 
large fluctuations, as we shall explain below. An even more problematic point is 
the tacit assumption that all moments of the Poincare recurrences are finite before 
taking the thermodynamic limit. The first explicit introduction of such a Poincare 
recurrence spectrum is done in jl7t . On the basis of numerical computations and 
heuristic arguments, it is claimed that the Poincare recurrence spectrum and the 
Renyi spectrum do not coincide even in the setting of Bowen-Gibbs measures : They 
argue that the Poincare recurrence spectrum must behave like 1/q when q — > —oo, 
which is not the case for the Renyi spectrum. Our goal is to clarify this claim in 
view of practical estimation of these spectra. 

We will mainly concentrate on hitting times because they are simpler to analyse 
and, at the same time, do share the same properties with return times for strongly 
mixing measures like Bowen-Gibbs measures. Our tools are thcrmodynamical for- 
malism and a very sharp result 1 that gives the error term in the convergence of 
P|) to the exponential law both in the size of An and in t. From this we derive two 

^The symbol b„ ^ c„ means that max(6„/c„, c„/fe,j) is bounded from above. We shall use it 
in the sequel. 
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approximation results: a global one and a local one. Indeed, Theorem 13.11 gives an 
approximation, for any n, of 



■wl^dfix^ forallqeR (8) 

as certain partition functions; Theorem I8. 21 gives an almost-sure approximation of 

\og{w.n{x,y)fi{An{x)) . (9) 

From (jSl we deduce our large deviation results and an explicit formula for the 
hitting-time spectrum, whereas from Q we derive a central limit theorem and even 
a law of iterated logarithm. 

As a matter of fact, we shall see that (1/n) logw„(a;, y) has the same normal 
fluctuations as —{l/n)logfi{An{x)), but their large deviations are not the same in 
some region. This is because 

wl dfix^^Yfi{Ary~'', forq>-l 

dnXfi^Y^{An)'^, forq<-l. 

This behavior is numerically observed in 117! for return times. 

Concerning return times, we can completely analyse lognormal fluctuations but 
not get an explicit formula for the return-time spectrum. Namely we prove that 
this spectrum coincides with the Renyi spectrum on [0, -l-oo). The reason for this 
is the presence of 'too soon recurrent cylinders'. Nevertheless, at the end of the 
paper we study non-overlapping return times f„. For them we prove that J d^ 
becomes flat for q < —1 and coincides with the Renyi spectrum for q S [0, cxd[. 
We conjecture that, in fact, the return-time spectrum really coincides with the 
hitting-time spectrum for Bowen-Gibbs measures. 



Outline of the paper. In Section [21 we record relevant definitions and results 
on hitting times as well as on Bowen-Gibbs measures. In Section |31 we establish 
the two main theorems of the paper for hitting times, namely a strong global ap- 
proximation of the 'free energy' of hitting times of n-cylinders, and a strong local 
approximation of hitting times. In Section^ we derive a number of corollaries from 
these two theorems: pointwise convergence, precise large deviation estimates, an 
explicit formula for the hitting time spectrum, a central limit theorem and a law of 
iterated logarithm. In Section|Slwe deal with return times. SectionElcontains three 
subsections. One is about (non-overlapping) return times. We improve our previous 
results by considering non-overlapping return times. Another one is concerned with 
bibliographical notes and possible straightforward extensions of our work. The last 
one illustrates that for the Manneville-Pomeau intermittent map, the hitting-time 
and return-time spectra are infinite for q > qc, where qc > (but we have a finite 
invariant measure). 



2 Set-up and background 

The phase space X will be the set of sequences x = {xi,X2, ■■■) where Xj E A (the 
finite alphabet), that is X = A^'^. The dynamics is given by the shift map T defined 
as {Tx)j = Xj+i for all j > 1. We only consider full shifts since the passage to 
subshifts of finite type is not an issue. Given a" 0102. ..a„, Oj G A, we denote by 
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[a"] the corresponding cylinder set that is [a"] = {x G X : xj — aj, j = 1, 2, n}. 
A point X X defines a sequence a cylinders that we naturally denote by [a;"], 
n > 1. The notation xf, 1 < i < j, stands for 2:^x^+1... Xj. The natural cr-algebra 
we take is the cr-algebra generated by cylinder sets. We omit to mention it in the 
sequel since it will always be the reference cr-algebra. 

Definition 2.1. We define the (first) hitting time of x to a cylinder [a"] as follows: 

V.](a;)^=*'inf{j>l:xf"-i=al'}, 

and the following hitting time, given x,y G X : 

w„(x,2/) t[,jj(2;) = inf{j > 1 : yf ""^ = .<} 

which is nothing hut the (first) time one sees the n first symbols of x appearing in 
y, i.e., the first time that the orbit of y hits the cylinder [x"]. 

The time w„(a;, y) is also called the waiting time |29| . 

Let us record the useful facts on the class of ergodic measures we are interested 
in. We refer the reader to (7| |^ 02 for full details. Let the potential Lp : X ^ R 
be of 'summable variations'. This means 

var„(^ < 00 

n>l 

where var„(^ sup{|(/?(x) — ^{y)\ ■ x'l — yi,x,y G X}. The condition imposed in 
[2] is more restrictive since it is var„(^ < (76*", for some C>O,O<0<1 (Holder 
continuity). 

BowEN-GiBBS PROPERTY. Assume that (fi : X ^ R has summable variations. 
Then there is a unique shift-invariant measure 11 = fi,p, that we call a Bowen-Gibbs 
measure, such that for all n > 1, for all a" and for any x G [a"] 

^ , < K (10) 

exp(-7ii^((p) + Snf[X)) 

where K = K{ip) > 0, Sn^{x) = ip{x) + (p{Tx) + ■■■ + ip{T''-^x) and P{ip) is the 
topological pressure of (p. From lfTU|l it is easy to deduce that 



1 

n — *oo 77, 



P(^) = lim^log^exp(S„^(a?)) 



where 5„(/?(a") sup{S'„(y9(x) : x" = a"}. We assume without loss of generality 
that (fi is normalized, which implies thatP(iy9) = and ip < [32] ■ A Bowen-Gibbs 
measure can also be characterized as an equilibrium state, that is the (unique) shift- 
invariant measure 77 that maximizes J dr] + h,j (T) , the maximum being equal to 
P{(p). This is the Variational Principle. Since we assume P{(p) = 0, this leads to 

- J ^ dfi = h^{T) . (11) 

For any g e M, define 

Miq)"^ lim -logV ^(K])i-« = lim - log / ^i{[x'(]y^ d^i{x) . 

n^oo Ji ^ n — !-oo Ji J 
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Using H1U|) . we trivially have that 



M{q) = lim - log / e~'?^"^(^' d^l{x) 



(12) 



n — >oo Ti 



It can be easily showed that q i-^ M.{q) is a well-defined function on M when is 
a Bowen-Gibbs measure. Moreover this function is convex and increasing. Indeed, 
by using the Bowen-Gibbs property (|10|l and the definition of topological pressure, 
we easily get M{q) = P((l - q)ip) (recall that P{(p) = 0). 

We now state the key result allowing us to analyse fluctuations of hitting times. 
The first result (with its proof) can be found in 0. 

Key-lemma 2.1 (Exponential distribution of hitting times with error 
term). Assume that ^ is a Bowen-Gibbs measure. Then there exist strictly positive 
constants c, C, pi, P2, with pi < P2, such that for all n £ N, all cylinder [a"] and all 
t > there exists p(a") G [pi,P2] such that one has 



Remark 1. The previous result is established under the hypothesis of -mixing of 
the process. Bowen-Gibbs measures indeed have this property. When the potential is 
Holder continuous, the proof of this property is done in '71 and in fact the ^-mixing 
coefficient decreases exponentially fast. When the potential has summable but not 
exponentially small variations, the ijj-mixing property is established implicitly in the 
proof of Theorem 3.2 in ^321. Notice that in this case we do not know how fast the 
ijj-mixing coefficient decreases. 

3 Strong approximations of hitting times 

We can state the main theorem of this section. For sequences of real numbers (6n), 
(c„), the notation 6„ ^ c„ precisely means that max{6„/c„, c„/6„} is bounded from 
above. 

We have the following 

Theorem 3.1. Let fi be a Bowen-Gibbs measure. Then 




(13) 



where £{a'^,t) = C exp(— cn) exp(— p(a")t). 




(14) 



and 




for g < — 1 . 



(15) 



Proof. Let q> Q. Then 




(16) 




(17) 
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By Key- lemma 12.11 there exist positive constants A, B such that for any t > Q one 
has 

Kev-lemma 12 . 1 1 also easily gives the lower bound : 

/ (t[,„] > dt>K'- Cexp(-cn) if" 

where < K' f'^ e-^^* < c» and < K" = f-^ e-P^* dt < oo. 

There exists an integer no such that for all n > hq, K' — Cexp(— crt) if" > 0. 
Therefore we get 



where 



Ki^q {K' ~ Cexp(-cno) isT") , K2 ^ qA f'-^e-^^dt. 

Jo 



This establishes (|14|l for q> 0. (The case (7 = is trivial.) 
Let now q e (—1, 0). 

J w-l^l dfi X M - Em(K]) / t-^l?l d/x 

= l9lE/^(K])^+l'l r i-l'l-V|tK] < TTT^I^^- (18) 
We first obtain a lower bound for the integral in the last expression : 

Hence, what matters is only the behavior for "large t" . Using again Key-lemma 
12.11 we get 

For all n > ni, where rii = i log + 1, we have 1 — (1 + Ce^™) e^''^ > 0. 

Since 1 - (1 + Ce-'^") e'P^* >\ - (1 + Ce-"") e'P^ for all i > 1, we obtain, for all 
n> ni, 

We now turn to the upper bound. We obviously have 

The integral from ^ to cxd is finite. Now we observe that, for every < i < ^, we 
have the following estimate: 

This estimate follows from the following lemma which is found in ^ Lemma 9]. 
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Lemma 3.1. For any integer t with t/i([a"]) < 1/2, one has 

l0g^i{t[al^] >t} 

Pi < , nn ^ P2 . 

where pi, p2 are the (strictly positive) constants appearing in Kev-lemma \2.1\ 

Using this estimate for the integral running from to i we get a finite upper 
bound since 

1 - e^P^*^ 



dt < oo . 



Hence we conclude that 



where K[ and K2 arc strictly positive constants. Hence we obtain (|14|) for q G 
(-1,0). 

Finally, let us consider the remaining case 9 < — 1. Then for sufficiently large n 
(such that ^([a"]) < 1/2), one has 



r /"^/^ 






/ + 








Jl/2- 






+ l2{ 


",<)]■ 



Clearly the second integral /2("-, a") is uniformly bounded in n. Indeed, 

.00 ^ 

However, the first integral /i(n,a") is diverging when n — > 00. Therefore the 
limiting behavior as n —s- 00 is determined by 

l'zlE/^(K])'+'"^i(^,<) = 

We again use Lemma [3. II to get 

provided that i < i. Using (HUJl (and P(<y9) = 0) we obtain 

R-^ eyip{-c'n) < pi[a'l]) < K e^p{-cn) 
where c, c' > 0. Hence we get 
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where we used the fact that for all k G R, 1 — e"*^ < k. Notice that for n large 
enough, the term between parentheses is strictly positive. Now, using the fact that 
1 — e~'^ > k/2 for any k e [0,1], and remembering that pi/2 < 1 (•^), and using 
again the Gibbs property (|10|l . we obtain 

where the term between parentheses is strictly positive provided that n is sufficiently 
large. Therefore, for n large enough, we end up with 

|g|pi(l-2l«l-iii:e-™) ^ / wn^'^^dfixfi ^ 2|g|p2(l-2l'?l-^/i:-^e-^") 



2(191-1) - EajMK])^ 

(Notice that L'Hopital's rule shows that there is no problem at q = —1.) Thus, we 
obtain (|15|) . which finishes the proof. □ 

We now turn to local strong approximation estimates. 

Theorem 3.2. Assume that ii is a Bowen-Gibbs measure. Then there exists Cq > 1 
such that for any e > eo , one has 

-elogn < log(w„(a;,y)/i([x'i'])) < loglog(n') (19) 

eventually p-a.s. (^). 

Proof. We want to find a summablc upper-bound to 

^ixfi{{x,y) : log(v^r„(a;,y)^(K])) > \ogt} = 

^A*(K]) /.{log(t[,.]/i(K])) > \ogt} (20) 

"■1 

where t will be suitably chosen. 

First observe that the function e(a", t) < Ce""^" for all t > in H13|l . Throughout 
this proof this error bound will be sufficient for our purposes. Using H13() in (|20|) 
yields 



HX^i{{x,y) : log(w„(a;,y)/i(K])) > logt} < Ce" 



e 



-pit 



Take t — tn — log(n'^), where e > is to be chosen later on, to get 

^ix^l{{x,y) : log(w„(a;,y)M([xi])) > loglog(n^)} < Ce"™ + . 

Choose e > l/pi. An application of the classical Borel-Cantelli lemma tells us that 

log(w„(x,?/)/i([a"])) < loglog(n'') eventually a.s. . 

For the lower bound, observe that using (|13|) with the same simplified error bound 
as before, we get for alH > 

fixp{{x,y) : log(w„(x,y)M(K])) < logt} < Ce'^" + 1 - e""^* < Ce'^" + ^2^- 

Choose t = t„ = n^"^, e > 1, to get, proceeding as before, 

log(w„(x, ?/)/i([a;"])) > — elogn eventually a.s. . 

The proof is finished by observing that both bounds hold for any e > max(l, pi). □ 



^Indeed, p2 = 2, see Uj- 

^By "eventually fix ^-a.s." we mean that there exists a set G with iJ,Xfi{G) = 1 and such that 
for any z £ G there is an integer TV = N{z) such that for all n > iV the inequality holds. 
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4 Corollaries 

In this section we derive the corollaries of Theorems 13.11 and 13.21 



4.1 Almost-sure convergence 

The following result tells us the 'typical' behavior of w„(a;,y). Throughout, h^{T) 
is the (Kolmogorov-Sinai) entropy of (X, T, /i). 

Corollary 4.1 (Almost-sure convergence of hitting times). Let ^ be a 
Bowen-Gibbs measure. Then 

hm -\ogWn{x, y) = hJT) ior fix ^ - a.e. {x,y) . 

n—^oo Tl 

Proof. By Theorem 13. 21 we get 

Um ilogw„(a;,?/) = hm -- log /x([x"]) = /i„(T) /xx^-a.e. 

n— >cx) Tl n^oo Tl 

where the second equality is given by Shannon-McMillan-Breiman Theorem. □ 

This result means that if we pick up randomly and independently of one another 
X and y, then the time needed for the orbit of y to hit [cc"] is typically of order 
exp{nhfi{T)). In fact, this result is valid under the more general assumption that 
the process is weak Bernoulli (or /3-mixing). For the details, we refer to |29j. Bowen- 
Gibbs measures are weak Bernoulli processes, see |7] for the proof. 



Remark. In j8| the author assumes that x is picked up randomly according to 
an ergodic measure r] whereas y is randomly (and independently) chosen according 
to a Bowen-Gibbs measure fi. The previous result becomes: 

lim - \ogWn{x, y) ^ h„{T) + hTivlfJ-) for?7X ^ - a.e. (a;, y) 

n— »oo Tl 

where hTirilfx) is the relative entropy of rj with respect to fi. The results obtained 
in this paper could be suitably generalized to that situation. 



4.2 Large deviations and multifractal spectra 

In this section, we study large deviations of ;i- logw„ around the entropy h^{T) 
where ^ is a Bowen-Gibbs measure, that is, we only assume that Lp has summable 
variations. 

The ansatz consisting in replacing /^([x"]) in the Renyi entropy by l/vfn(x,y) 
leads to the following definition. 

W{q) = lim -log f ^l{x,y) d{fix ^l){x,y) , (21) 

n— ►oo n J 

provided the limit exists. We do not use exactly the same definitions as in TT*. The 
present definitions are motivated by large deviation theory. 
Introduce, for convenience, the following functions: 

W„(g) = ^log J vfl{x,y) d{^ixfj,){x,y) 

for all n > 1 and g £ M. (Notice that Wniq) can be infinite.) 
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Observe that each function q i-^ yVnil) is a convex (increasing) function on R 
(hence, in particular, a continuous one 
Observe also that 

There is no Kacformula for hitting times for general ergodic dynamical systems. 
Ergodicity only ensures that almost surely there is a finite first hitting-time. The 
only fact we know without any assumption is that J t^^„j dfj, < oo for all q < 0. 
Indeed, 

for any g < 0. From [H] it follows that for any Bowen-Gibbs measure and for any 
cylinder [a"], we have 

d/i < oo for all g G M . 

Hence yV„(g) < oo for aU g G K, n > 1. 

We now turn to large deviation results. We refer the reader to ^I] for back- 
ground on this topic. 

Corollary 4.2 (Scaled generating function of hitting times). Assume 
that fi is a Bowen-Gibbs measure. Then 

fM(g), for g > -1, 

"'<"=ip(2;), J,-<-,: <^^' 

If fj, is not the measure of maximal entropy, then the function q s- yV'(g) is strictly 
convex on (— l,oo). 

Proof Clearly ^ and ijTHl imply □ 

Notice that g i~+ yV(g) is continuous (as it must be for a convex function on M) 
but not differentiable at g — —1. Indeed, it can be easily checked that the right 
derivative at —1 of W is not zero but equal to ~ J (p d^2Lp > 0, where fi2<p is the 
Bowen-Gibbs measure for the potential 2ip. 



Corollary 4.3 (Large deviations of w„). Let ^ be a Bowen-Gibbs measure 
which is not the measure of maximal entropy. Then for all u > we have 

lim -\og{nx n) \ -log^,, > h^{T)+u\ = inf {-{h^{T) + u)q + W{q)} 

and for all u G (0, uo), = | limqj^_i yV"(g) — h^j_{T)\, 

lim -log(^x^)(llogw„</i^(r)-?/l = inf {-(/i (T) - ?/)g + >V(g)} 

Notice that uq > hi^i{T), that is, we capture the large fluctuations of logw„/n 
above and below hi^i_{T) since yV'(O) = h^{T) (see the appendix for the proof). 



^Remind that a convex function defined on a finite interval of the real line can be discontinuous 
only at the endpoints of that interval. 
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Proof. By Theorem 13.11 and (|12|) we immediately get that for any q > —1 
W{q) = M{q) ^ lim - log / exp{-qSnf ) dfi . 

n^oo n J 

It can be easily deduced from |3J that the function q i-^ P{q'f) is and strictly 
convex if (and only if) /i is not the measure of maximal entropy. In the Holder 
continuous case, it is real analytic and also strictly convex if (and only if) ^ is not 
the measure of maximal entropy pl| . 

We can apply a large deviation result due to Plachky and Steinebach |2E1- (Recall 
that a strictly convex differentiable function has a strictly increasing derivative.) □ 

Let us remark that when the measure is the one of maximum entropy, there are 
no large fluctuations which is not surprising. 

The Renyi spectrum is defined here as M(g) =^ Jv[{q)/q and the hitting-time 
spectrum W{q) = W{q)/q. We get h^(T) for the value of these spectra at q — 
(using L'Hopital's rule). 

Corollary 4.4. For any Bowen-Gibbs measure we have the following: 

W(q)=M(g) forg>-l and W{q) = P{2ip)/q forg<-l. 



4.3 Log-normal fluctuations 

The purpose of this section is to show that Wn{x,y) and l//i([a;"]) have the same 
lognormal fluctuations for Bowen-Gibbs measures associated to Holder continuous 
potentials. Namely, we prove a central limit theorem and a law of iterated logarithm. 

We refer the reader to for full details on the central limit theorem for Holder 
continuous observables with respect to Bowen-Gibbs measures with a Holder con- 
tinuous potential. Define the following variance: 

cr'^ == lim — / ( Sn^ — n [ ip dji] dfi . (23) 
n^oo n J \ J J 

It is well-known that ii > (u^ < oo because of the exponential decay of 
correlations) and 

VteM lim mI^^^^^^^^^^^ <4 -A^(0,l)((-oo,i]) (24) 

where A/'(0, 1)((— oo,t]) — exp(— ^^/2) df. This is just a particular instance 

of the central limit theorem for Bowen-Gibbs measures where the observable 
is -Lp. By (dJ J -If dfi = hf,{T). One has ct^ = if and only if -p + h^{T) (or 
equivalently ip — J ip d/i^) is a coboundary, i.e. a function of the form Q — goT, for 
some measurable function g. This means that cr^ = if and only if /x is the (unique) 
measure of maximal entropy. 

Corollary 4.5. Assume that fi is a Bowen-Gibbs measure with a Holder contin- 
uous potential which is not the measure of maximal entropy. Then 

lin, ^ X ^ ( -nh,{T) ^ 1 ^ ^^^^ l)((-oo, ^]) . (25) 

Moreover, 

= Mm - [ (logw„ - hf,{T) f d{^ixfj.) . (26) 

n^oo n J 
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Proof. Our goal is to show that the central limit theorem for logw„(x,y) results 
from the one for — log/x([a;"]) which in turn results from the one for ~Sn^{x). The 
latter assertion is trivial since, for all a;, 

-Sn^ix) -C< -log^(K]) < -Sn^ix) + C (27) 

where C > 0, due to Bowen-Gibbs inequality H10() . Hence the quantities — log /i([a;"]) 
and —Snf(x) have the same mean h^i{T) and variance ct^. 

We now use the strong approximation formula (|19|l from Theorem 13 . 21 together 
with inequalities H27|l to get 

^ -r= — {) tor fix fi— almost every [x,y) . (28) 

By a basic result of Probability Theory (see ^21 for instance), the preceding result 
and H24|l imply the desired statement. 

The proof of l|26|l is given in the Appendix. □ 

We emphasize that (|26|) does not follow from Theorem 13. 21 see the Appendix. 
We now state and prove a law of iterated logarithm for logw„: 

Corollary 4.6. Assume that n is a Bowen-Gibbs measure with a Holder contin- 
uous potential which is not the measure of maximal entropy. Then 

logw„ - n/i^(r) 

limsup — 1 fixfi — a.e. (29) 

n->c» o■^/2?^loglogr^ 

Remark that we get —1 instead of 1 when taking 'liminf instead of 'limsup'. In 
fact, we could show that the set of accumulation points of the sequence {(logw„ — 
nhf^{T))/^/2n log log n}„ is the interval [— ct, +0"]. 

Proof. Using Theorem 13.21 we get that eventually jix /i-almost surely 

£ log n ^ logWn - nh^{T) ^ -Snf - nh^{T) ^ log log n" 
a\J2n log log n ~ a\/2n\o% logn a^/2n log log n ~ a\/2n\og \ogn 

Taking the limit supremum n — > oo and using the law of iterated logarithm for 
—Sn^, we finish the proof. □ 

The law of iterated logarithm for —Sn^p can be found in |12j . 

5 Return times 

We now turn to return times. As we shall see, we obtain less complete results than 
for hitting times. 

5.1 Set-up 

Definition 5.1. The (first) return time of a point x into its n-cylinder [x^], n>l 
is defined as: 

r„(x) = inf{fc > 2 : = Xi} . 

The following result is proved in 2, Section 6]. In order to state it we need to 
define the set of n-cylinders with 'internal periodicity' p < n: 

Spin) = {K] : min{fc G {1, ...,n} : K] n T-^[a^] ^ fb} = p} . 

Notice that the set of n-cylinders can be written as the union Ui<p<n '^p("')- 
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Key-lemma 5.1 (Exponential distribution of basic return times). Let fi 
be a Bowen-Gibbs measure. Then there exist strictly positive constants Ci,C2,C3 
such that for any n G N, any p € {1, ■■■,n}, any cylinder [a"] € Sp(n), one has for 
all t > p 



^r'''"^^^^ CKRKI) ' K]|-C«)exp(-0 <£(ar,i), (30) 

where e{ai,t) = Ci exp(— C'2n) exp(— Csi) and^(a") is such that] C(oi)— p(ai)| < 
D exp(— C27i); with D > and p(a") given in Kev-lemma \2. 1\ Moreover, 

fi{z : r„(z) > t I [a"]} 1 for all t < p. 
5.2 Large deviations 

We define the following (possibly infinite) quantities, for all e M, provided the limit 
exists: 

Jl{q)'^ lim ilog / rliz) d^i{z) (31) 



^n{q) = - log / r«(z) d^(z) and 



Without any assumption on /i, the function Ji trivially exists at g = and equal 
zero. If fi is assumed to be ergodic then Dl(l) — log \A\. Indeed, 

Cr^{z) dn{z) ^ ^ (^j Vn{z) (z)^ M(K])- 

By Kac's formula the integral is equal to l//x([a"]), hence we get 

lim —log / r„(z) d^(z) = hm — log^^ 1 = log |^|. 

n — >oo 77, y n — ^oo 77, ^ — ^ 

a? 

Let us emphasize that the finiteness of 31 for all g > 1 is not obvious at all. We 
need to know the finiteness of Jj^„j d^ for all g > 1, that is the finiteness of the 
moments of return times to [a"]. This point seems to have been overlooked before. 
For q < 1, all moments of return times are of course finite due to Kac's formula 
(T-invariance is indeed sufficient) but nothing rules out a priori the possibility that 
the moment of the return time be infinite beyond a certain order q > \ for some 
TiQ (and hence for all n> since the moment of order g > as a function of n is 
increasing). This will be illustrated in Section IHTI From [5] it follows that for any 
Bowen-Gibbs measure and for any cylinder [a"], we have 



/ d/^[a"] < oo for all q G 



Hence %i{q) < oo for aU q e M, n e N. 

Proposition 5.1 (Partial large deviations for r„). Let ^ be a Bowen-Gibbs 
measure which is not the measure of maximal entropy. Then 

lim ilog^i j-logr„ >h^{T)+u\ ^ mi {-{h^,{T) + u)q + ']l{q)} 
n— ►oo n yn J (;>0 

that holds for all u > 0. 
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Proof. From Kev-Lemma 15 . ll we can deduce that 

/'r^dM^^MK])^"^ (32) 

for all g > (the case q = is trivial). This implies that M(g) — 3?(<z) for all g > 0. 
Thus we can again apply Plachky-Steinebach's large deviation estimate on K+. We 
leave the details of the proof of H32|l to the reader since it is very similar to the one 
for hitting times we gave above. □ 

Some large deviation results are given in lO for r„, but they are valid only in 
a small (non-explicit) interval around hf^{T). This restriction is due to 'too soon 
recurrent cylinders'. Because of the same problem, we can only extend this result 
to the whole range of possible deviations above h^{T). 

5.3 Lognormal fluctuations 

We summarize what happens for return times in the following theorem and corollary. 

Theorem 5.1. Assume that fi is a Bowen-Gibbs measure. Then there exists ei > 1 
such that for any e > ei, one has 

-elogn < log(r„(x)M(K])) < loglog(n^) (33) 

eventually fi-a.s. 

Corollary 5.1. Let fi be Bowen-Gibbs measure associated to Holder continuous 
potential. Then Corollaries \4.5§4-(>] hold for r„ instead o/w„. (/ix/x is replaced by 
H.) 

Sketch of proofs. The proof of the corollary follows exactly the same lines 
as that for hitting times. Let us sketch the upper bound, leaving the lower bound 
to the reader. To apply the classical Borel-CantcUi lemma, as before, we need to 
upper bound 

M{x:log(r„(x)/i(K]))>logt I [<]} = 

n 

53 53 MK])M{log(r„MK]))>logi I K]} (34) 
where t will be suitably chosen. We now use Kev-Lemma 15 . II to get 

n 

m<J2 E ^(Wl) ((P2 + D)e-'^' + Cie-<5-) 

P=l [a'i']e5p(n) 

for alH > p /x([a"]), where we used the fact that there is some uq such that for all 
n > no, pi/2 < pi- < (^(^^n^ < p2 + D. Now choose t = t„ = \og{n'), 

e > and notice that i„ > p /^([a"]), for all p = 1, ...,n for n large enough since 
m([«i]) ^ K by the Bowen-Gibbs property. Therefore we get 



which is summable in n provided that e > 2//3i. We leave the rest of the proof to 
the reader. 
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6 Final comments and open questions 



6.1 Hitting-time and return-time spectra in presence of in- 
termittency 

We consider here the following intermittent map (the so-called Manneville-Pomeau 
map) defined on the interval [0, 1]: 

T : X x + modi, a £ (0,1) . 

The thermodynamic formalism for such a map is now well-understood. We refer 
the reader to the recent paper |22| for more details and references on what we will 
use. 

Let if = — log|T'|, the potential function. The map T admits an absolutely 
continuous invariant measure fj, which is an equilibrium state for the potential if. 
(fi is not the only equilibrium state. Any measure of the form tii + {I — t)6Q, where 
t E [0, 1] is an equilibrium state for (p; Sq is the Dirac measure at 0.) This means 
that P{ip) — hfj {T) + J dji, where P{'^) is the topological pressure of ip and 
h^{T) is the measure-theoretic entropy. In fact P{ip) = because of the Rokhlin 
formula. It is well-known that M((7) = P{{1 — 0)1^). From the behavior of the 
pressure function we get the following properties for JA^q): it is continuous, convex 
and non-decreasing. Moreover, P((l — q)'p>) = for g < 0, P((l — q)<f) > for g > 
and q ^ -P((l — q)^p) is a real-analytic function for q > 0. At the critical point one 
has the following asymptotics: 

^«i^^Mr)as,iO. 

q 

The Manneville-Pomeau map has two intervals of monotonicity, /o, Ii, from 
which one can define cylinder sets: Iii,i2,...,i„{x) = H T~^Ii^ n • • • n T^^Ii^ is 
that interval of monotonicity for T" which contains x. Contrarily to the case when 
maps are everywhere expanding, the ratio 

^(-^il ,22, in (-'')) ii^) 



is not uniformly bounded in n and x. This comes from the fact that distortions 
are not bounded. A more careful analysis shows that one can find bounds from 
above and below which are polynomial in n and uniform in x. Such a measure is 
an example of a weak Gibbs measure. 

The following basic proposition shows that large deviation results in the usual 
sense do not hold. 

Proposition 6.1. For all q > ^, we have for all n > I that Wn{q) = 00 (hence 
yV{q) = 00, but M(g) is finite for every g S K. The same occurs for 3?n((z) for all 
q> - + 

Proof. We are going to show that 'Wi{q) ~ log J wf d{iix /i) becomes infinite from 
some go = 9o(q^) > on (hence yVn(g) = 00 for all n when q > qo since Wn+i^x, y) > 
Wn{x,y) for all {x,y).) The main point is the following estimate for /x(/oo...o): 



where C is a positive constant. But /i(/o...o) (with i symbols) is nothing but the 
measure of points that do not enter the right interval Ji before £ iterations. Now 
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observe that 



— £ symbols 

Therefore Wnio) = oo for all n> 1 and q> 

To pass from Wn{q) to 3^„((7) use Proposition 1 in □ 

It was recently proved in that, for < a < ^, 1/ ^{li^^i^ i^i^)) r„(a;) 

have the same lognormal fluctuations 

Proposition Vi . II shows that the return and hitting time spectra arc not relevant 
for non-uniformly hyperbolic systems. Indeed, a single indifferent fixed point makes 
these spectra infinite for all q > qda). Moreover, it implies that 1/ ^J■{Iil,i2,...,ir,i^)) 
and and r„(a;) cannot have the same large deviations. 



6.2 More on return times 

In |17j the authors study the recurrence-time spectrum ^{q). They show heuris- 
tically that 3?(g) must behave like 1/q as q oo. (Be careful of the different 
convention used therein to define 3?(q).) Two numerical simulations confirm this 
behavior: The graphs of the spectrum really look like a constant divided by q for 
g > 1. This is indeed what we get rigorously (remember that q = 2 corresponds 
to g = —1 in our convention) for hitting-time spectrum W{q). We are not able, 
as we saw above, to prove this for return times. The difficulty comes from 'too 
soon recurrent' cylinders. More technically speaking, we do not have the analog 
of Lemma 13.11 for return times and therefore the finiteness of the integral in (|18() 
escapes us. 

We can go a bit further by looking at non-overlapping return times f „ (studied 
recently in 0). By definition, such return times cannot be 'too small': The (first) 
non-overlapping return time of a point x into its n-cylinder [x"], n > 1 is defined 
as: 

f„(^) =inf{fc>l:xittf =xn. 
We have the following approximation result: 

Proposition 6.2. Let ^ he a Bowen-Gihhs measure. Then, for every q < —1, 
For every q > 0, 

(The symbol ^ is precisely defined at the beginning of Scction|3|) Therefore, if 
we let Ji{q) be the analog of ^{q) where f„ replaces r„, the previous result implies 
that 

■3i{q) = P{2(p) foraUq<-l and ■k{q) ^ M{q) for aU q > 
(Use the Bowen-Gibbs property and the definition of topological pressure of Section 

Recall that for hitting times we proved a more precise result fTheorem l3.1f) since 
the second approximation works for q S [— l,oo[ in that case. 
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Proof. First write 

/oo 
rld^l = Y,k-\^\ ^l{r,,^k}. (36) 
fc=i 

Using the i/j-mixing property (see for details) we get that 

m(K])(1 - i'ik - 1)) < MK]{f« ^k}< m(K])(1 + ^(fc - 1)) . 

Hence, 



2 

Therefore we get 



1^ ^ /f?. ^3 



where Cq = X^feLi ''^ < 0° since |(?| > 1. We leave the proof of the statement in 
the range q G [0, oo[ to the reader (use the analog of Kev-lemma l5.1l for f„ that can 
be extracted from 0). 

So, we arrived at the desired result. □ 

We conjecture that Theorem 13.11 is true with return-times instead of hitting- 
times. It could be easier to first prove this conjecture for non-overlapping return 
times. 

6.3 Relevance of the hitting-time and return-time spectra 

In view of |14j (saturation of level sets) and what we obtained in the present paper, 
one can legitimately ask what is the relevance of the return time spectrum, except 
as a trick to compute the Renyi spectrum for g > —1. The same can be said for the 
hitting-time spectrum. Even in the comfortable setting of Bowen-Gibbs measures, 
these spectra contain no information for q < —1. In presence of intermittency, we 
saw that they are infinite for q > qc- 



6.4 Related works and an extension 

It is worth to indicate to the reader the differences between our work and the pre- 
vious ones. In the paper |1U| . the authors study only return times for Bowen-Gibbs 
measures associated to a Holder continuous potential. They prove a central limit 
theorem and a large deviation principle. Here we not only improve the lognormal 
approximation but also extend the range of accessible large deviations above the 
true entropy. Moreover we handle potentials with summable variations and not only 
Holder continuous ones for large deviations. We also mention PU] for a general (but 
much less sharp) result on lognormal fluctuations for return times. In the context 
of i/j-mixing stochastic processes, there are two references [201 IS]. Both deal with 
local strong approximations, in particular central limit theorems. The author of |2(J| 
directly uses the ■i/'-mixing property. The author of |34j first proves an approxima- 
tion to the exponential law of rescaled hitting and return times and then deduces 
strong local approximations. We emphasize that his approximation is much less 
sharp than the one we use here. This difference is not relevant for deriving strong 
local approximations but becomes essential to handle large deviations. The only 
paper dealing with large deviations of hitting times is P] where the authors study 
the first occurrence of a cylindrical pattern in the realisation of a Gibbsian random 
field on the lattice Z'^, d> 2. Our proof is very similar to that of this work. 
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We also note that our results can be extended to a more general class of processes, 
namely the processes satisfying the (^-mixing property with a summable (/s-mixing 
sequence. This is because Key- lemma 12 . II is proved not only for ^/;-mixing but also 
for such processes 1 . But up to our knowledge, this does not define a natural class 
of equilibrium states on shift spaces. That is why we did not state our results under 
this assumption. Nevertheless, an interesting class of non-Markov expanding maps 
of the interval has this property with an exponentially decreasing (/j-mixing sequence 
(with respect to the partition given by the discontinuity points of the map). This 
class was indeed studied in 23_, . We could therefore sharpen the results of that paper 
since Corollarv l5.ll apply. We could of course write down the analogous results for 
hitting times. Regarding large deviations of hitting times, we could derive some 
approximations in the spirit of Theorem 13 . II and derive some estimates like that of 
CoroUarv 14.31 But one has to be careful with the control of some 'bad' cylinders 
for which the "distorsion property" (the analog of Bowen-Gibbs property (|10|l in 
that context) does not hold, which is the price to be paid for the non-Markovian 
partition. 



7 Appendix 

We prove the convergence in mean and in quadratic mean of logw„. The former is 
related to the slope of W at 0, see the comment after Corollarv 14.31 The latter is 
related to the proof of H26|l in Corollarv l4.5l We emphasize that Theorem 13 . 21 does 
not help because almost sure convergence does not say anything for convergence 
unless we have bounded random variables, which is definitively not the case here. 

It is easy to get that yV,'i(0) = n/^^S^n dfix ^. Since q i-^ Wn{q) is a 
convex and, at least, continuously difierentiable function on (— l,+oo), yV"(0) — 
lim„^oo >V4(0). Let us show that hm„^oo W;(0) = h^{T). 

First observe that convergence in mean of (logw„)„ is equivalent to showing 

lim - I \\og{w„{x,y)fi{[x'^])\ dnxfi{x,y) = 0. 

n-+oo n J 

Indeed, 

^ J \ \ogw„{x,y) -nh^{T)\ d^x^i{x,y) < 

|log(w„(x,y)/i([x5'])| d^ixfi{x,y) + ^ J \ log fiHx'^]) + nh^{T)\ d^i{x) . 

The second term goes to zero by using ()10|l. P{(p) = and (|11|1 and noting that 
/ Sn^ d^ ^ n J ip d^ hy T-invariance of /i. We have 

log(w„(a;,2/)M([<])| dfixn{x,y) V^([a'i']) / ^{log(t[„.]/i(K])) > t'} dt' 
The change of variable log t = t' leads to 

/°° df 
A.{t[,,.]>VMK])}y 

By Key-Lemma 12.11 this integral is finite and bounded between, say, Ci and C2 
(independent of n). Therefore we get the desired result. 

Now turn to the convergence in quadratic mean of (log w„/r7,)„. Observe that 
the following identity holds: 



^ J {logwn ~ nhf,{T)f d{fiXfi) 
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- / (log(w„MK]))) rf(MXA') +- / {log y.{[xl]+nh^{T)fdii 



log(w„M(K]))(logM(K]) d{^,x^,) 



The second term goes to (see formula (|23ll and use pU|) '). Hence the proof is 
done if we show that the two other terms go to as n — > cx). Proceeding as before 
we get that the integral in the first term equals 

2 EA^(K])/°°M{tK] > VMK])} ^ dt. 

Using again Kev-lemma l2.1l we bound the integral from above and below uniformly 
in n. Now consider the integral in the third term which is equal to 

J2 KK]) i^ogKK])+ nh^{T)) /'log(t[,.](2/V(K])) dt,{y). 



We recognize the same integral as above which we know bounded from above and 
below uniformly in n. The factor in front of the integral was also treated above in 
this section. The proof is finished. 

We could prove the same results for r„ by using Key- lemma 15.11 and for f „ by 
using the corresponding result found in ,5,. 
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